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properties 


Lal Babu Sah Telee™, Ram Suresh Yadav’, Vijay Kumar 


ABSTRACT 


Here, a new distribution having three called Modified Inverse Lomax Distribution is 
proposed. Important statistical properties like the survival, hazard rate, the 
probability density function (PDF) is studied. Least Square, Cramer-Von Mises 
and Maximum Likelihood estimation methods are used for estimation of 
parameters using R programming software. A data set is discussed and 
performed the goodness-of-fit to assess the application of the proposed 
distribution. Various methods of model comparison and model validation are 
also used. The proposed model called Modified Inverse Lomax Distribution is more 
applicable as compared to some existing probability model. 


Keywords: Lomax distribution, Estimation, Hazard function, Cramer-von Mises, 


maximum likelihood. 


1. INTRODUCTION 


Probability distribution helps to simulate the real-life problems and to compute 
the real-life data precisely and effectively. During last decade, different new 
probability distributions are introduced by researchers using different techniques. 
These techniques may me adding some extra parameters to distribution, merging 
the distribution or inverting the variables etc. These methods make new 
distribution more flexible and useful than the existing distributions. Many 
standard techniques to analyze the problem of real-life data are available; we still 
need models to solve the problem more effectively and precisely. That is in all the 
cases, classical techniques are not effective as the new distributions in literature, 
we can find numerous distributions. We can derive parametric distribution by 
changing the number of parameters to existing distribution such as Lomax and 
exponential family of distribution (Marshall and Olkin, 1997). Marshall and Olkin 
family of distribution Ghitany et al., (2007) is used to extend Lomax distribution. 
We can find different new distributions derived by using the Lomax 
distribution. Power Lomax distribution is more flexible than existing Lomax 
distribution with decreasing and inverted bathtub hazard rate functions (Rady et 
al., 2016). Taking alpha as power, a new distribution named alpha power inverted 
exponential distribution was introduced using inverted exponential distribution 
(Ceren et al., 2018). Using Lomax random variable as a generator Ogunsanya et 
al., (2019), Odd Lomax exponential (type III) distribution was introduced. 


Compounding of odd generalize exponential using inverted Lomax distribution 
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for obtaining Odd generalized exponentiated Inverse Lomax distribution (Maxwell et al., 2019). Lomax exponential distribution is 
obtained using Lomax distribution (Ijaz and Asim, 2019). The inverse Lomax- exponentiated G- family Falgore and Doguwa, (2020) 
is introduced using Inverse Lomax distribution. In real life, we can find many data that have bathtub-shaped hrf. In literature we 
can also find many modifications of Weibull distribution. Two parameter Weibull distributions is given as 


F (x, a, 8) = exp[-(a, x)’ 


This distribution is modified to generate several distributions that possess bathtub hrf. Modifications of Weibull distribution is 
exponentiated Weibull distribution (Mudhokar and Srivastav, 1993). As given by Lai et al., (2003) to get new lifetime distributions 


as 
F(x) = exp[ax’ .exp(Ax)] 
We have modified the Modified Lomax distribution to introduce new probability model called, Modified Inverse Lomax 


distribution where inverse of the variable is taken. CDF and PDF of two parameters Lomax distribution (Pathak and Chaturvedi, 
2013) is, 


F(x)=1-(1+ Bx)” 3x>0,(B,4)>0, 
oa f(x, B.A) =AB(1+ Bx)? 5x>0,8>0,4>0 


An extra parameter a is added to modify the Lomax distribution as Modified Inverse Lomax distribution. CDF and PDF proposed 


model is given as. 
F(x) =[1+ (Bie ]", x>0, a,B,4>0. 
And the PDF of Modified Inverse Lomax Distribution is given as 
f(a =a(B/ xe) tAne* [14 (Blane |; x>0, a8, >0 


The Modified Inverse Lomax (MILX) distribution 
A three parameters Modified Inverse Lomax distribution with CDF and PDF are as follows 


Fay=|i+(B her , x>0, a,B,A>0. (1) 


Xx 
And pdf of MILX is, 


f(x) = a( 4) [1 Bie je" i (l+Axje“*; x>0, a,B,A>0 (2) 
x 


Survival function 
Survival function of MILX is 


S(x)=1-[14(B/ae* |" x>0, a, B.A >0 ) 


Hazard rate and reversed hazard rate function 


MILX model has hrf as 
n(x)= FO), 0<x<0 
1= F(a) 
=a(B ix [1 +(Biae"] d+eane|1-[14+ Blane" | @ 


Reverse hazard function of MILX is 

f(x) 

F(x) 
=alpleyleane™[Lepive*]'[eGineT 6 


hey (x) = 
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The various shapes of pdf and hazard rate function of MILX (a - B . A) for at various values of parameters are displayed below 


in (Figure 1). 
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Figure 1 Graphs of PDF (left) and Hazard rate function (right) for some values of parameters 


Here we have discussed the quantile function also. Quantile function can be defined here by expression (6) where X is non negative 


random variable. 
QO, (p)= | (p) 
log B -logx—Ax+1- pp” =0 (6) 


Skewness and Kurtosis 
Quartile based coefficient of skewness defined by Bowley’s is, 


_ Q(3/4)-2Q(1/2)+Q(1/ 4) 
a Q(3/4)-Qd/4) 
Also, Moors, (1988), defined coefficient of kurtosis based on octiles can be defined as 
Q(0.875)-Q(0.625)-Q(0.125)+Q(0.375) 
Q(0.75)-Q(0.25) 


K(Moors) = 


2. PARAMETER ESTIMATION 


Maximum Likelihood Estimation 
Here, we have presented the ML estimators (MLE's) of the MILX model are estimated by using MLE method. Let a sample 


x= (x, pees X, ) drawn randomly from MILX (a 7 B ; A) with log likelihood function, 


£ = nloga +nlog B+ > log (1+ Ax,)-2 log x, -A>x, 
i=l i=l 


i=l 


-(@ +1) log[1+(B ix), ee] 
i=l 
Differentiating equation (7) with respect to parameters as 
a ot — Dog (1+(B/ xje* ) 


Ca a 4 


(7) 
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AX; 


aL x a 
eeu < 


op =I X; a 


of n n exp(— AX) 
— = B(at+l)+ 
A Dis: TFA 5 DT (Bix,)e** | 
Equating a = ie ad =0 and solving simultaneously for all parameters, ML estimators of the MILX(@, B, A) 
da OB OA 


distribution can be obtained. Normally, it is not possible of solving non-linear equations above so with the aid of suitable computer 


software one can solve them easily. Let © =(a@ ,B,A) denote vector of parameters of model MILX (@ ‘ BA) and the 
a Beh -1 
corresponding MLE of © as@ =(@,f,A) thus asymptotic normality becomes, (0 -@) > N, |0.(7 (O)) | I (OQ) is 


the Fisher’s information matrix defined as, 


(= | al | al 
6a? Caop Cada 
1( = al (=) al 

Opoa op OB OA 


2 2 2 
gi so Vf et.) alot 
adda) | adap an? 


-1 
Generally, 0 is not known. That is, MLE has an asymptotic variance (I (9)) is worthless. Asymptotic variance can be 


approximated by plugging in the estimated value of parameters. We have used O(Q) as an estimate of [ (0) and is given by 
om al al 
0G = AaOB AaOA 
2 2 2 
Gey 2S) sae 
apoad af? afar 
a’l a’l a’l 
ACG =OAOB OA? 


laa) 
Matrix H is called Hessian matrix. 
Newton Raphson algorithm is used to optimize Likelihood function. Matrix below is the called Variance Covariance matrix, 


4 var(@)  cov(a, B) cov(a. , r) 
[1 (9), =| cov(B, 4) var(B) —cov(B, A) (8) 
cov(A, a) cov(A,B) — var(A) 
Let Z,,, is upper percentile of SNV then using asymptotic normality of MLEs, approximated 100(1-b) % CI of @, B, and A of 


MILX(a@, 8,4) is defined as, 
+ Z, 19 / Var(@) , B+Z,,,\vat(B) and A + Zap var(A) . 


Least Square Estimation 
We have also used the LSE to estimate the @, B, & A of MILX. For this we have to minimize function below 
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A(x|@, B.A) = ¥] FW) ~ | (9) 


i=l 


Let us consider that F(X,,)) is CDF of order statistics x < Mig Sve Bey taking 1s >, Career X,} as a sample having n 


(1) 
items from a distribution function F (.). LSE of parameters are calculated by minimizing the function 


A(x|a@, B,A) ay {(£}e"} — (10) 
i=l (i) 


Differentiating (10), we get, 


CA gt a | | B) m=) B\ = \\- 
bao 1 eae (i) _ l 1 abe, (i) 1 aa (i) 
Hef} aan) (ele) 
= ~(a+l) 

A pis 7 (1 exp(—Ax,;,) Bb 7 
op ad eas | a) | Xi) Wel pas 

-~a —(a+l) 
La ne Lay ore se Ei weer 
a = 2aift Y exp AX) eas in| waa eas ix 


We can obtain weighted least square estimators by ee 


B(X;a,B,A) ->W, FX, jae 


n+1] 


1 (n+1) (n+2) 
Weights Wi are W, = ——— = +—_~~___ 
© Var(X)) — i(n-i +1) 


Thus, weighted LSE of a, 8, and A are given by minimizing the function below, 
2 


B(X;a, B.A) = 5 EY (1+ (Braue yr -— (tt) 


‘a i(n-i+l) n+l 


Cramer-Von Mises estimation 


This method of estimation can be obtained by minimizing the function 


Z(X:a, B.A = Lad] [o+(ara em yr" a 


42 

1 a | = ‘} 2i-1 

=—+ 1-(1+ Bxe — 12) 
12n | ( Px ) 2n | 


Differentiating (12), we get, 


< = 25 {( + (2 Xj )e Je Aree yr - = 1] log(1 ~ (2 thes je ) (1 + (2 t Tis eon er 


aZ ning WO" _ 26-11 ~-y “ 
F --2a}|{(1+(87 0 )e ) 7 i [< Ge ()) 
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-a-l 


OZ a ite, we Sie BAI ae 

rie wo lf4 (Bray )e*} - —] ie | 
OZ a 
da OB OA 


After solving non-linear equations = 0 CVM estimators can be obtained. 


3. APPLICATION TO REAL DATASET 


For applicability of the model, we have considered a real set data. Data is by Lee and Wang, (2003) (Statistical methods for survival 
data analysis, 3rd edn. Wiley, New York). Ramos et al., (2015) (The Kumaraswamy-G Poisson family of distributions. Journal of 
Statistical Theory and Applications, 14(3), 222-239). 

0.08, 2.09, 3.48, 4.87, 6.94, 8.66, 13.11, 23.63, 0.2, 2.23, 0.26, 0.31, 0.73, 0.52, 4.98, 6.97, 9.02, 13.29, 0.4, 2.26, 3.57, 5.06, 7.09, 11.98, 
4.51, 2.07, 0.22, 13.8, 25.74, 0.5, 2.46, 3.64, 5.09, 7.26, 9.47, 14.24, 19.13, 6.54, 3.36, 0.82, 0.51, 2.54, 3.7, 5.17, .28, 9.74, 14.76, 26.31, 0.81, 
1.76, 8.53, 6.93, 0.62, 3.82, 5.32, 7.32, 10.06, 14.77, 32.15, 2.64, 3.88, 5.32, 3.25, 12.03, 8.65, 0.39, 10.34, 14.83, 34.26, 0.9, 2.69, 4.18, 5.34, 
7.59, 10.66, 4.5, 20.28, 12.63, 0.96, 36.66, 1.05, 2.69, 4.23, 5.41, 7.62, 10.75, 16.62, 43.01, 6.25, 2.02, 22.69, 0.19, 2.75, 4.26, 5.41, 7.63, 17.12, 
46.12, 1.26, 2.83,4.33, 8.37, 3.36, 5.49, 0.66, 11.25, 17.14, 79.05, 1.35, 2.87, 5.62, 7.87, 11.64, 17.36, 12.02, 6.76, 0.4, 3.02, 4.34, 5.71, 7.93, 
11.79, 18.1, 1.46, 4.4, 5.85, 2.02, 12.07 

By use of R software of the optim () function Team, (2020), calculation of MLEs of MILX by maximizing the log likelihood 
function defined in equation (5.3.1) (Mailund, 2017). Values determined are tabulated in (Table 1). Figure 2 displays the graph of 
profile log-likelihood function of parameters showing that ML estimates is uniquely defined. 


Table 1 Estimated parameters using MLE, LSE and CVME 


Method | alpha | beta lambda 
MLE 1.0701 | 6.7270 0.0589 
LSE 0.8818 | 10.8745 | 0.0882 
CVME_ | 0.9001 | 10.6555 | 0.0899 
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Figure 2 Profile Log-likelihood function of parameters. 
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Graph of P-P plot and Q-Q plots are shown in (Figure 3). It is found that MILX model has better fitting to the data taken in 
consideration (Figure 4). MLE, LSE and CVE methods are used for estimation and the estimated value of the parameters of MILX 
model with log-likelihood, AIC, BIC, CAIC, and HQIC criteria in (Table 2). KS, W and A2 statistic and calculated p-value are given 


in (Table 3). 


Theoretical distribution function 


0.0 


0.2 


04 


06 0.8 


Empirical distribution function 


Figure 3 The P-P plot (left) & Q-Q plot (right) of the MILX Model. 


60 


40 
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20 


Table 2 Estimated parameters with LL and values of information criteria 


Empirical quantiles 


Method | alpha | beta lambda | LL AIC BIC CAIC HQIC 
MLE 1.0701 | 6.7270 | 0.0589 | -401.084 | 808.1674 | 816.7235 | 808.3609 | 811.6438 
LSE 0.8818 | 10.8745 | 0.0882 | -402.620 | 811.2404 | 819.7965 | 811.4340 | 814.7168 
CVME_ | 0.9001 | 10.6555 | 0.0899 | -402.759 | 811.5186 | 820.0746 | 811.7121 | 814.9949 
Table 3 The KS, A2 and W statistic and their p-values 

Method | KS(p-Values) | A2(p-Values) | W (p -Values) 

MLE 0.0433 (0.9702) | 0.3678(0.8797) | 0.0535(0.8556) 

LSE 0.0347(0.9978) | 0.3538(0.8928) | 0.0297(0.9776) 

CVE 0.0346(0.9979) | 0.3774(0.8706) | 0.0291(0.9793) 
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Fitted quantiles 


x Sample quantiles 


Figure 4 Histogram and density function of fitted distributions (left) & fitted quantile vs sample quantile (right) of different 
estimation methods of MILX. 


In this part of study, applicability of MILX taking a real dataset used previously is presented. Model’s potentiality is compared 
with following four distributions. 
I. Generalized Exponential (GE) distribution: 
The Pdf of GE model (Gupta and Kundu, 2007) is 


for (x) = ale {I : gaye 


II. Exponentiated Exponential Poisson distribution (EEP): 
Pdf of EEP is, (Risti¢ and Nadarajah, 2014) 


apa _ ae Pe 
Cee ’*(—e**)*exp{-A(l-e ae 
(l=e™) 

III. Generalized Exponential Extension (GEE) distribution: 

Pdf of GEE is (Lemonte, 2013) 


Sore =a BA (1+ Ax) soe [1 expfi—(+aay"}) exp{I-(I+ 2.x) “| 


IV. Lindley-Exponential (LE) distribution: 
The PDF of LE (Bhati et al., 2015) distribution is, 
2 


Fie y= 4 = {1-log(1—-exp(-Ax))} exp(-Ax)(1-exp(-Ax))°"; 


Various information criteria values for the testing of the applicability of the MILX are tabulated below (Table 4). We have 
displayed the graph of goodness-of-fit of MILX and models defined as (Figure 5), We have plotted the empirical distribution 
function of the proposed model MILX and the fitted distribution function and is shown in (Figure 6) 


Table 4 AIC, CAIC, BIC, HQIC and log likelihood Values 

Dist. AIC. CAIC BIC HQIC ll 

MILX | 808.1674 | 808.3609 | 816.7235 | 811.6438 | -401.0837 
LoW | 824.9487 | 825.1423 | 833.5048 | 828.4251 | -409.4744 
EEP 825.5056 | 825.6991 | 834.0617 | 828.9819 | -409.7528 
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GEE | 827.2026 | 827.3961 | 835.7586 | 830.6789 | -410.6013 
LE 828.0985 | 828.1945 | 833.8026 | 830.4161 | -412.0493 
GE 830.1552 | 830.2512 | 835.8592 | 832.4728 | -413.0776 
1.0 
— Empirical 
0.6 
a 
Oo 
04 
02 
0.0 4—— 
T T T 


F(x) 


Figure 6 Empirical distribution functions versus the fitted distribution function of the MILX 


80 


Goodness-of-fit of the MILX model with other competing model are compared in this section. We have also displayed value of 
Kolmogorov-Simnorov test, Cramer-Von Mises and, Anderson-Darling test statistic in (Table 5). MILX has the minimum test 


statistic values. That is, p-values of MILX model are higher indicating the better and consistent fit than considered models. 


Table 5 Test statistics and p-values 


Dist. | K-S(p- Values) | W (p - Values) | A2 (p - Values) 
MILX | 0.0433(0.9702) 0.0535(0.8556) | 0.3678(0.8797) 
LoW_ | 0.0321(0.9994) 0.0149(0.9997) | 0.1010(0.9999) 
EEP 0.0380(0.9925) 0.0220(0.9946) | 0.1486(0.9987) 
GEE 0.0441(0.9640) 0.0393(0.9370) | 0.2631(0.9630) 
LE 0.0622(0.7059) 0.0899(0.6377) | 0.5250(0.7211) 
GE 0.0725(0.5115) 0.1279(0.4652) | 0.7137(0.5472) 
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4. CONCLUSION 


Presented article explains a distribution named Modified Inverse Lomax distribution having three parameters. Detailed analysis of 
different statistical characteristics of model is also explained. We have presented expressions for its hazard rate function, survival 
function, the quantile function and skewness & kurtosis. Maximum likelihood estimation, Cramer-Von Mises estimation, and Least 
Square estimation methods are used to estimate the parameter. It is found that MLEs are better w.r.t. LSE and CVM. PDF curve of 
model has shown that it can have various shapes like increasing as well as decreasing. Monotonically increasing, constant as well as 
bathtub shaped based hazard function is seen. Applicability and suitability of model is evaluated. For this, we have considered a 
real-life dataset. We found distribution that MILX is much flexible. 
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